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Abstract 



We study waves in a rod of finite length with a viscoelastic constitutive equation of fractional 
distributed-order type for the special choice of weight functions. Prescribing boundary con- 
ditions on displacement, we obtain case corresponding to stress relaxation. In solving system 
of differential and integro-differential equations we use the Laplace transformation in the time 
Oh domain. 
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1 Introduction 

ON , Fractional derivatives have been used in describing physical phenomena such as viscoelasticity, 

diffusion and wave phenomena. There are two approaches in formulating differential equations 
with fractional derivatives in physics and mechanics. In the first approach classical "integer order" 
differential equations of a process are modified by introducing fractional derivatives instead of 
integer order ones (see books by Mainardi 1997, Podlubny 1999, and Kilbas et al. 2006). In the 
second approach one uses variational principles such as the Hamilton principle as a starting point for 
deriving equations of a process, where a modification of the classical case is achieved by replacing 
some (or all) integer order derivatives in Lagrangian density by fractional derivatives of certain 
kind. Then the resulting Euler-Lagrange equations are equations of a process and they contain both 
left and right fractional derivatives (see papers by Agraval 2002, Atanackovic & Stankovic 2007, 
Atanackovic et al. 2008). 

In this paper we generalize classical wave equation for one-dimensional elastic body by following 
the first approach. Recall the classical setting. Consider the equation of motion 

—a(x,t)=p—u(x,t), jce[0,L], t>0, (1.1) 
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where p, a and u denote density, stress and displacement of a material at a point positioned at x 
and at a time t, respectively. It is coupled with the Hooke Law 

a (x, t) = EE (x, t) , 16 [0,L], t>0, (1.2) 

where E is a modulus of elasticity and £ is a strain measure, defined by 

d 
£(x,t) = - 7 -u(x,i), x£ [0,L], t> 0. (1.3) 

ox 

Combining (|l.l|l - (|1.3|) , the classical wave equation is obtained as 

d 2 p d 2 

—u(x,t) = ——u(x,t), x£[0,L],t>0. 

We propose the generalization of a constitutive equation (|1.2|) by replacing it with a constitutive 
equation which corresponds to a generalized viscoelastic body: 

fa (a) Q D?a(x,t) da = E <f> 2 (a) D?£ (x, t) da, ie[0,L],t>0, (1.4) 

Jo 

where E is a positive constant (having dimension of stress), <j) l and 2 are given functions or distri- 
butions and oD"y is the left Riemann-Liouville fractional derivative of a function y £ AC ([0, T]) , 
for every T > 0, of the order a £ [0,1), defined as 

where T is the Euler gamma function. Recall, AC ([0, T]) denotes the space of absolutely continuous 
functions (for a detailed account on fractional calculus see a book Samko et al. 1993). In case 
when <p l and <p 2 are distributions, we assume that <p 1 and (f> 2 are compactly supported by [0, 1] 
(0i, 02 ^ £' (K) , supp 0i, supp 2 C [0, 1]). In this case integrals in (jl.41) are defined as 



(a) D?h (t) da, ^ (t) ) := (0 (a) , ( D?h (t) , ^ (t))) , <p £ V (R) . 

supp / 

For details see a work by Atanackovic et a/. (2009a). Recall, T>' + (R) denotes the space of distri- 
butions supported by [0, oo) and (h (t) , ip (i)) denotes the action of a distribution h £ T>' + (R) on a 
test function ip £ T> (R) (see a book by Vladimirov 1984). 

In (|1.4j) . 0! and 2 denote constitutive functions or distributions that are determined experi- 
mentally. The constitutive equations of type (|1.4j) were used earlier in papers by Hartley & Lorenzo 
(2003), Atanackovic (2002o), Atanackovic et al. (20096) and Atanackovic et al. (2009c). There are 
number of forms that 0i and 2 can take (see a paper by Hartley & Lorenzo 2003, for example). 
In the sequel we assume that 

0i(a):=a Q , 2 (a) := b a , a £ (0, 1) , a < b. (1.5) 

The restriction a < b follows from the Second Law of Thermodynamics (see for example papers by 
Atanackovic 2002a and Atanackovic 2003). If a = b, then (jl.41) reduces to the Hooke Law. The 
choice of 0i and 2 in the form (| 1 . 5|) is the simplest choice guaranteeing dimensional homogeneity. 



Note that with <j) x (p) := 8 (p) + rfS {ll — a) and (j) 2 {ll) := E^t^S {ll — /?) (6 denotes the Dirac 
distribution) we obtain 

o + r%D?a = EooT^Dfe, (1.6) 

while with <f) l (ji) := 6 (ll) + t*S (jj, - a) and <f> 2 (/x) := E (S {ll) + t£<5 (/z - a) + t%5 {ll- /?) ) we 
obtain 

<r + r^i}^ = So f 1 + r«o^r + ^o^f) e. (1.7) 



Recall, system (|1.1|) . (|1.3|l and (jl.6|l . respectively system (jl.ip . (II .3|) and (jTTTJ) , was treated in a 
work by Rossikhin & Shitikova (20016), respectively in a work by Rossikhin & Shitikova (2001a). 
Also note that the distributed order dissipation of type (jl.41) was also used in the context of one 
degree of freedom mechanical systems in papers by Atanackovic et al. (2005) and Atanackovic & 
Pilipovic (2005). 

Our aim is to find functions u and cr, locally integrable on R and equal to zero for t < 0, so that 
these functions satisfy (|1.1|) , (11.31) . (|1.4|l . for x G [0, L] and £ > 0, as well as the appropriate initial 
and boundary conditions. Actually, we will introduce dimensionless quantities and transform the 
system ([i~Tj) . (fl~3| . (fi~4|) , subject to (|2~Tj) and (|2T2|) . into the system (|2~3|) . subject to (|2~4|) and 
(I2~5|) . 

The paper is organized as follows. In SJ5] we introduce dimensionless quantities, proceed by 
formal calculation and by the use of the Laplace transformation we obtain solutions to (|1.1[) . (|1.3[) . 
(|1.4p in the convolution form. We impose initial conditions as well as boundary conditions to 
(jl.ip . (jl.31) . (|1.4p . Boundary conditions describe a rod that is fixed at one of its ends, while 
the other end is subject to a prescribed displacement T (this is the case of stress relaxation if 
T = To-ff , with H being the Heaviside function) . Section [3] is devoted to the calculation of the 
inverse Laplace transformation, which leads to the explicit form of a solution. More precisely, we 
investigate some properties of functions in order to be able to apply the Cauchy residues theorem, 
which is used to calculate the inverse Laplace transformation. We obtain displacement u and 
stress a for the boundary condition T = T H in §3.1.11 as well as for T = T n H + F, where F 
is an appropriate function supported by [0,oo) , in §3 .1.21 We conclude that solutions are locally 
integrable functions supported by [0, oo) . Moreover, they are smooth functions for t > 0. Numerical 
examples corresponding to stress relaxation are presented in 2J Concluding remarks are given in 



2 Convolution form of solutions. Formal calculation. 

We prescribe initial conditions for system (jl.ip . (jl.31) . (jl.4p 

d 
u(x,0)=0, —u (x, 0) = 0, a(x,0) = 0, £(x,0)=0, x£[0,L}. (2.1) 

We subject system (jl.ip . (|1.3[) . (jl.4p to boundary conditions corresponding to case of stress relax- 
ation 

w(0,t) = 0, u(L,t) = T(t), tet. (2.2) 

Function T is locally integrable function equal to zero for t < 0. 
Introducing dimensionless quantities 

_X-t_U_fJ-T (h- <^~ 



V LyfE ' " V- E 1 - V ^ ( L y|) Q ' V2 (Lyf)"' 



and using the fact that the fractional derivative transforms as 

Dfu(i) = (iy|j oD?u(t), 
we obtain, after omitting bar over dimensionless quantities, the following system 

^ (*>') = ^ (*'*)> 

(t) 1 (a) Dfa(x,t)da= (f> 2 (a) D?£ (x,t) da, (2.3) 

Jo 

£ (x, t) = — u (x, i) , x € [0, 1] , i > 0. 
ax 

System (|2.3p is subject to initial 

d 

u(x,0)=0, — w(x,0)=0, cr(x,0) = 0, £{x,0)=0, x G [0, 1] , (2.4) 

and boundary conditions 

u(0,t) = 0, w(l,*) = T(i), tel. (2.5) 

Additional assumptions on T will be given in JJ3.1.1I and JJ3.1.2I 

In the sequel we assume that (p-^ and <f> 2 are functions given by ()1 .5[) . Using (jl.5p in (|2.3p -? and 
formally applying the Laplace transformation to (|2.3|) and (|2.4[) . we obtain 

— ct(x, s) = s 2 u(x, s), 

a(x,s) / (as) a da — £ (x, s) \ (bs) a da, (2-6) 

Jo Jo 

d 

£ (x, s) = — u (x, s) , x e [0, 1] , s € D. 

Recall, the Laplace transformation of / <G L^ oc (K) , / = in (— oo,0] and \f (t)\ < ce at , t > 0, for 
some a > 0, is defined by 



/>oo 

/ (s) = L [f (t)} (s) := / / (£) e^di, Re s > a 
Jo 



and analytically continued into the appropriate domain D. Domain D for (|2.6p is determined after 
(j2~7l) . bellow. 

System (|2.6[) reduces to 

d 2 , , o In (6s) as — 1 . . _ 

— ^u x, s - s 2 — i— i- -u x, s = 0, x € 0, 1 , s € £>, 

ax^ In (as) bs — 1 

whose formal solution is 

u (x, s) = Ci (s) e xs V 7 ^!^ + C 2 (s) e " xs V^T tfiF a; e [0, 1] , s e L>, (2.7) 



where C\ and C2 are arbitrary functions which will be determined from the boundary conditions. 
Since the natural logarithm has the branch point at s — 0, we have D = C\ (— 00, 0] and this will 
be used in proposition 13. II 

Applying (|2.5[> i we obtain C = C\ = —C2 and thus 



u (x, s) = C (s) (e'VSi tt^ _ e ^VS^t j f x e [0, 1] , s e C\ (-00, 0] . 

Using (|2.5I) 9 in the previous expression, it follows 

u(x,s)=f(s)P(x,s), x G [0, 1] , s G C\ (-00, 0] . (2.. 

We introduced P as 



sinh (xs / ln ( 6s ) as - 1 
P(x,s) = - \ V^f 6 *- 1 / , zG[0,l], S GC\(-oo,0]. (2.<)i 

Note that P (x, •) is a distribution supported by [0, 00) . It is clear that for x = 1, P (l,t) = S (t) , 
i€t. Actually we will calculate P and show that for every x G [0, 1] , P (x, •) is a locally integrable 
function on K, equal to zero for t < 0. 

Since T and P are supported by [0, 00) , displacement u is given by 

u(x,t) = T(t)*P(x,t), x G [0,1], t£R, so that (2.10) 

u(x,t) = 0, i€ [0,1], t < 0, 

where we use * to denote the convolution. Recall if /, g £ L] oc (R) , supp f,gC [0, 00) , then 
(/ * 3) (*) : — Jn f ( T ) 5 (t ~ T ) dr, i € R. Calculation of (j2.10[) will be done by the use of the Laplace 
inversion formula applied to (|2.8j) . Moreover, we will show that uh (x, t) — H (t) * P (x, t) , t € R, 
is a continuous function, equal to zero for t < 0. 
Next, we use ()2.6j) 9 and (|2.6|U and obtain 



a (s, *) = ^(^j^_i ^" (*. s ) - s G [0, 1] , s G C\ (-00, 0] . (2.11) 

In order to determine <r, we use ()2.8j) and ()2.9[) in p. lip , so that 

cr (x, s) = sf (s) f (x, s) , a; G [0, 1] , s G C\ (-00, 0] , (2.12) 

where 



ln(bs) qg— 1 



-. , /in (as) 6s - 1 C0Sh ( XS \J \n(as) bs-1 , 

Ta: ' s = VT~k( 1 / /w / / -' ^e 0,1 ,SGC\ -oo,0. (2.13) 

Note that T is a locally integrable function for t > 0. For the determination of a, we will use the 
Laplace inversion formula applied to (|2.12[) and obtain 

a(x,t) = ^(T(t)*T(x,t)), are [0,1], £ G R, (2.14) 

where the derivative is understood in the sense of distributions. Again, a {x, t) = for x G [0,1], t < 
0. For the detailed account see 



3 Explicit forms of solutions 

In this section we will calculate inverse Laplace transformations of distributions and functions on 
R supported by [0, oo) . For that purposes, let us define 



"<H^i^' secu -~' 01 - (3 - i) 

In the sequel we will write A(x) ~ B (x) if lim -^l — 1. Next proposition establishes some 
properties of M. 

Proposition 3.1. (i) M is an analytic function in s € C\ (— oo, 0] ; 

(ii) lim M (s) = 1 and lim M (s) = 

s— s-0 |g|— s-oo 

sec\(-oo,o] sec\(-oo,o] 

(raj £e£ p € (0, sq) , so > 0. Then 



M(p±LR)~ ./j— 1— (/(In (oiZ) In (6i?)) 2 + f^ln-') e Ti arctan W^i 
V o in (ait J y \2 a J 

as R — > oo. 

Proof. We first prove (i). The only points where M could be singular are s = - and s = t. Since, 

ln M = ln ( l + (6 g -l)) = "(-!)" ( to -1)" _ 

6s -1 6s - 1 ^ n + 1 ' ' iv, J, 

n— 

it is obvious that s = ^ is a regular point of M. Similar arguments hold for s — -. Limits in (ii) 
can easily be calculated. In order to prove (iii), let us introduce /i = y/p 2 + R? and v = arctan—. 
It is obvious that /i ~ R, v ~ ^, as R — > oo. Then M (p ± ii?) becomes 

M (p ± ii?) 



'In (an) In (6/u) + ^ 2 T i^ln£ /(ap - 1) (bp - 1) + a6i? 2 ± LR (6 - a) 
In 2 (a//) + i/2 y (bp - l) 2 + (6i?) 2 

[a6i? 2 + (ap - 1) (bp - 1)] [In (an) In (6 M ) + v 2 \ +Rv(b- a) In £ 



((6i?) 2 + (6p-l) 2 ) (ln 2 (a M ) + ^ 



R(b-a) [In (a/i) In (b/j,) + v 2 } - v In \ [abR 2 + (ap-l) (bp - 1)] 



((6i?) 2 + (6p-l) 2 ) (ln 2 M + ^ 2 ) 



--— -—-Win (ai?) In (6i?)=Fi7r ln -> as i? ^ oo. 
b In (ait) V 2 a 



D 



3.1 Determination of displacement u in case of stress relaxation 

We investigate properties of P, given by (|2.9|) . and find a solution to (|2.3j) . ()2.4|) . (|2.5j) . with weight 
functions given by (|1.5j) . in two steps. First, we find a solution with the boundary condition ()2.5[) 
in the case when 

T (t) = T H (t) , T o >0,t€R. (3.2) 

Then we proceed to a more general case, when T is assumed to be of the form 

T (t) = T H (t) + F (t) , ten, (3.3) 

where F is a locally integrable function, equal to zero on (— oo, 0] . Additional assumptions on F 
will be postulated in JJ3.1.21 

Let us examine properties of P, given by (|2.9[) . Clearly, it has complex conjugated poles at 
pSn , n € N. Poles are solutions of 

sinh (sM (a)) = 0, i.e. sM (s) = ±nm. (3.4) 

Let us examine the position and the multiplicity of solutions to (|3.4p . 

Proposition 3.2. There are infinitely many solutions psh, , n € N, o/ (5^1), such that 



Re ( P5 (±)) « 4 " V ° y (3.5) 

In I yabnir J In b\/-nir I 

Im LW) « ±i? « ±J-nn, (3.6) 

as n — > oo. Moreover, there exists uq £ N, such that poles pSn , for n > no, are simple. 

Proof. Let us square (|3.4p and put psi ' = Re 1 ' 1 ', £ (— it, it) . Then, after separation of real and 
imaginary parts, we obtain 

R 2 cos (20) Re (M 2 (i?e i0 )) - R 2 sin (20) Im (M 2 (Re^)) = -nV, (3.7) 

i? 2 sin (20) Re (M 2 (i?e^) ) + R 2 cos (20) Im (M 2 (ite^) ) = 0. (3.8) 

By the use of p.ip , real and imaginary parts of M 2 (i?e'^) are 

(In (aR) In (bR) + 2 ) (a&i? 2 - (a + o) i? cos + l) 



Re(M 2 (i?e I, P)) = 



+ 



Im (M 2 (Re^)) = - 



(In 2 (aR) + 2 ) (6 2 i? 2 - 2W?cos0 + 1) 
ln^(6-a)i?0sin0 
(In 2 (aR) + 2 ) (6 2 i? 2 - 2W?cos + 1) ' 

01n£ (abR 2 - (a + b) i?cos0 + l) 
(In 2 (aR) + 2 ) (6 2 i? 2 - 2W?cos 0+1) 
+ (b-a)R sin (in (aR) In (6i?) + 2 ) 
(In 2 (aR) + 2 ) (b 2 R 2 - 2bRcos + 1) ' 



Letting R — > oo, previous expressions are written as 

abR 2 In (aR) In (6i?) a In (W?) 



Re(M 2 (i?e^)) 
Im(M 2 (.Re i0 )) 



b 2 R 2 In 2 (ai?) 6 In (ai?)' 



ao In ^i? 2 t; 



a & 
= — — In ■ 



1 



b 2 R 2 In 2 (aR) ~ b a^ln 2 (aR) ' 
Using (23), flU} and (13~TU1) . we obtain 

Im(M 2 (i?e^)) 



lni 

O 

Re (M 2 (i?e^)) ~ ^ln (oi?) In (bR) ' 



(3.9) 
(3.10) 

(3.11) 



tan (2(f)) 

Let € (0, tt) . Then ^^ > and 

Since <fi ^ and tan (20) — ► 0, it follows that <fi —> ^ from the interval G (?, ^) . Therefore, by 

(J3.11I) . we have 



tan(2</ 



-)• 



as i? ^ oo. Hence, </> € (0, | ) or cf> e (f , ^) 



sin d> fa 1 



Sin* 

2 a 



2 In (ai?) In (W?) / 
Inserting (JHSJ), (|3~T0f and (|3~T2|) in ([57f]l . we obtain 

1 



1 , COS ( 



^ln^ 

2 a 



21n(aR)1n(bR)' 



In 2 (aR) In 2 (&£)+(§ In £)' 



" ,2 " >fl ) + iM) 



i? 2 



(3.12) 



(3.13) 



Thus, real and imaginary parts of ps„ , as R — ► oo, obtained by (|3.12[) and (J3.13I) . are as stated 
in proposition. 

In order to prove that solutions to (|3.4p are simple for n > no, we define 



Then (s€ C\(-oo,0]) 
d 



f(s):=sinh(sM(s)), seC\(-oo,0] 
ln| (b-a)s 



d J(s) M (s) I 1 - 2 ^ ^ 1q ^ , 2 ^^ _ ^ ^ _ 1 1 



cosh (sM (s)). 



Solutions to f (s) — are given by (|3.4I) . and so, as 



psk 



— > oo, 



ds 



/(*) 



„(±) 



(-iy 



M(s) 1 



In- 



(& — a) s 



2 In (as) In (6s) 2 (as - 1) (bs - 1) 



c<±) 



By proposition 13.11 M ~ y^ and this implies that 

d 



els 



/(*) 



(-1)" 



p.S 



(±) 



oo. 



Thus, for large 



psi. 



we have 4-f (s)\ _ <±> 7^ and solutions are simple for n > n . 



D 



3.1.1 Case T = T n P 



This is the case that has physical importance, since we obtain displacement in case of stress relax- 
ation test. Formally, we write (J2.10J) with (|3.2[) as 

u H (x,t) = y H(t)*P(x,t), xe [0,1], i€R. (3.14) 

The following theorem is on existence and properties of uh- 

Theorem 3.3. Let T = Yo-ff and let 4> x and <p 2 be given by lil.5\) . Then the solution to 
\24\) , Ii2.5\) is given by {3.1$ , where 



1 f 00 / sinh (xqM (qc- i7r )) sinh (xqM (ge i7r )) 
2riJ 1 sinh (qM (qe~ i7T )) sinh (qM (ge i7r )) 



+ ^[Res(p(x, S )e st , P 4+) 

n=l 

+ Res (P (x, s)e st , P s ( n - ) 
P(x,t) = 0, xe [0, l], t< 0. 

The residues are given by 



- qt dq 



x € [0, 1] , t > 0, 



Res 



(PMe^psW) 



sinh (xsM (s)) 
A[ sin h( s M( S ))]' 



(3.15) 
(3.16) 



(3.17) 



.(±) 



and simple poles ps n , /or n > no, are solutions of \3.J$ . Function P is real-valued, locally 
integrable on R and smooth for t > 0. 
T7ie explicit form of solution is 

T f°° fsinh (xqM (qe- [n )) sinh (xgM (ge iw )) \ 1 - e -9 * 
UffOM) - — ^ y S i n h ( q M(qe-^)) sinh (qM (qe™)) j q ? 



/ f^[Res(p(x,s)e 



ST PS (+) 



|dr, x G [0,1], * > 0, 



+ Res(p(a;,s)e ST ,ps(r ) 
u H (x,t) = 0, ie [0, 1], i < 0. 
Function ur is continuous at t = 0. 

Proof. We calculate P (x,t) , x € [0, 1] , £ G R, by the integration over a suitable contour. 
Let £ > 0. The Cauchy residues theorem yields 

r. OO 

i P (jc, s) e s *ds = 27ri J^ fees (P (x, s) e st , P s ( r +A + Res (p (x, s) e s \ P s ( - 



(3.18) 
(3.19) 



(3.20) 



Im s 




Res 



Figure 3.1: Integration contour T 

where r = Ti U T2 U I"3 U T e U T4 U T^ U Tq U 7 , so that all poles lie inside the contour T (see hgure 
151) . 

First we show that the series of residues in Q3.15P is convergent. By proposition ^. 21 poles pSh, 
of P, given by (|2.9[) , are simple for n > n . Then residues in (I3.20[) can be calculated as it is given 
in ([5TT7|l . We use (JO) to write (|3~T7)l as 

Res(p(a;,s)e s *,p4 ±) ) 

sin (wrx) >•<■•' 



("1) T 



1 



ln- 



s(b-q) 



2 ln(os) ln(6s) r 2(as-l)(bs-l) _ 



n > riQ. 



„c±) 



Let ps; = Re ±1( ^ 7 then (|3.21j) transforms into 



(3.21) 



Res 



e st P s {±) 



(P(X )S ) 

„ sin (mrx) 



(-1)' 



Dpi?i cos 0p±i(0+i?i sin 0) 



I ITT 



hi 



s(6— a) 



21n(as)ln(6s) ' 2(as-l)(fes-l) 



5 = ^0*'' 



= (-1) 

and therefore, for n > uq, we have 
Res 



„ sin (n-KX) Re Rt cos * [cos (0 + Rt sin 0) ± i sin (<p + Rt sin <j>)] 

' mr l\ In 4 , sffe-ol 



21n(os)ln(6s) ' 2(as-l)(6s-l) 



S^floi'' 



(p (x, s) e st , P «W) + Res (> (*, s) e st , P 4r ) 
„ sin(n7ra:) m C os<p 



cos (d> + Rt sin <i) + i sin (<i + Pi sin < 



1 



ln- 



+ 



s(fe— a) 



2 ln(as) ln(6s) r 2(as-l)(6s-l) 



s=flc ir 



cos (6 + Rt sin 0) — i sin (d> + Rt sin ( 



1- 



lu 



s(b— a) 



21n(as)ln(bs) ' 2(as-l)(6s-l) 



s=Rc-"i> , 



(3.22) 



10 



Let n — > oo (then also 



P s 



(±) 



— > oo, i.e. R — > oo). Then 



1 



In- 



+ 



s (b — a) 



2 In (as) In (6s) 2 (as - 1) (6s - 1) 
and (|3.22p . as n — > oo, becomes 



=i£e± 1 * 



->1, 



Res (p (a:, s) e st , P jW) + Res (p (x, s) e st , P s^ 



sin (nirx) 



Re Rt cos ^ cos U + Rt sim 



tin 



Proposition 13.21 (|3.5p . implies that 
Re ( PS W) , 



7T 6/6 

— In — \ —n 
4 a V a 



lnfV^6n7r")ln[6y / f 



< -CVn, n > 



n . 



Also, by (|3.6p . we have that -^ ps \/ f 7r - This implies that summands in (|3.20j) can be estimated by 

Ke~ Ct ^™, which implies the convergence of the sum of residues in p.20p . 

Second, we calculate the integral over V in p.20p . Consider the integral along contour IV Then 



P(x,s)e st ds 



< 



P(x,p + iR) 



,( P +iR)t 



dp. 



Let R — ¥ oo. In order to estimate 



P (x,p ± iR) , using (iii) of proposition 13. 1[ we write 



M(p±\R) ~ v±iw, 

fa 1 In (aR) In (bR) 



b In (aR) 
a 1 



^J (\n (aR) In (bR)) 



2 + (fln|) 2 



fin* 

z a 



Then, as R 



P(x,p±iR) 



b ln (ai?) y(ln(ai?)ln(6i?)) 2 + (fln^)' 



sinh [x (pv — Rw) ± ix (pw + -Rv)] 



< 



sinh [(pv — Rw) ± i (piy + ito)] 

~x(pv — Rw) _i_ „ — a?(pv— _Ruj) 
g— (pv— -Rw) I 

I + c - 2x (P v ~ Rw ) 

I I — c —2(pv — Rw) I 



„-(l-a:)(pu-Hiu) 



0. 



(3.23) 
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The previous statement is valid since, as R — > oo, 

fa 1 1 



pv — Rw 



V bln ( aR ) ^(\ n (aR)\n(bR)) 2 + (l\n±y 

x (pin (aR) In (bR) + R- In - 
\ 2 a 



r± J^m+^L 



R 



, by]j In {aR) 2 a hi (ai?) ^/ln (ai?) In (W?) ^ 
Therefore, according to (|3.23[) . we have 



— > oo. 



lim 

R— ¥00 



P{x,s)e st ds 



= 0. 



By the use of p.23[) . we conclude that similar arguments are valid for the integral along the contour 
Tq. Thus, 

lim / P{x,s)e st ds = 0. 
r-hx J T6 

Next, we consider the integral along contour T2 

2xfle i *M(flc i '*) _ 1 



P {x, s) e st ds 



r-2 



< R 



1 R(l-a;)e i *M(fle 1 *) 



e 2_Re i *M(fl.c i *) _ 1 



^Rt cos0 



Since M ~ y^ as \s\ — > 00 and cos <j> < for <j) € [^, 7r] , by the Lebesgue theorem, we have 

P{x,s)e st ds 



lim 

R— >oo 



r 2 



< lim / i?e fico8 ^(* +(1 - a:) ^)d^=()^ 

R-Kx> 



Similar arguments are valid for the integral along the contour T 5 . Thus, 



lim 

R— Kx> 



P (jc, s) e"*ds 



r r , 



The integration along contour T £ gives 



lim 

£->0 



P(a;,s)e s *ds 



= lim 



={l-x)e i *M{e& i *) 



_ -2Kec i *Af (ec'*) 



1-e 



et cos r 



^ _ c -2ec i <t>M(ec i *') 

Integrals along parts of contour T 3 , T 4 and 7 give 

sinh(xgAf (ge i7r )) 



= 0. 



lim / P{x,s)e st ds 

R— foo /p 



sinh(gM(oe i7r )) 

f ~ , . f°° sinh (xqM (qe~™)) 

lim / P{x,s)e st ds = - / , ,\ „, ■ // 

^-ir 4 io smh(gM(ge— )) 

lim / P{x,s)e st ds = 2niP{x,t). 



e- qt dq, 
e' gt dq, 
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Now, by the Cauchy residues theorem (|3. 201) . the function P is determined by (|3.15p . 

In order to see that P is a real-valued function, we use M (qe ±17T ) , q € [0, oo) , and note that 
M (qe~ 17r ) = M (qe 17r ), where the bar denotes the complex conjugation. Due to the exponential in 
the hyperbolic sine, we have sinh (xqM (qe~ 17r )) = sinh (xqM (qe 17r )) and therefore the integrand 
in (|3.15j) is of the form 

sinh (xqM (qe'™)) sinh (xqM (qe™)) 
sinh (qM (qe~' m )) sinh (qM (qe i7T )) 

_ f sinh (xqM (qc™))\ sinh (xqM (qe i7r )) 



\ sinh (qM(qe i7T )) J sinh (qM (qe [% )) 
( sinh (xqM (qc i7T ))\ 



-2ilm 



I sinh(qM(ge i7r )) 



which implies that the first term in (13.15[) is real. 

Next, we examine Res ( P (x, s) e st , pSn J in order to prove that the sum of residues is also real. 

By (ET2TT) and 



In- 



(b - a) s 



2 In (as) In (bs) 2 (as - 1) (bs - 1) 



1 



In- 



(b — a) s 



2 In (as) In (bs) 2 (as - 1) (bs - 1) 
we obtain that Res (P (x, s) e st , ps n j = Res (P (x, s) e st , ps„ 



,(+) 



It is clear that 



Res (P (x, s) e st , P s^ +) ) + Res (P (x, s) e^^s^A = 2 Re ^Res (P (x, s) e st , P s 



( + ) 



This implies that the second term in (|3.15[) is also real for n £ N. Hence, (|3. 15|) is a real- valued 
function. 

Let t < 0. We prove that the integral over 7 does not depend on the choice of So (see figure 
13. ip . Let f = 7 U 7j U 7q U 7 2 (see figure |3~^|) . where so and s' are chosen so that all poles, i.e. 
solutions of (13.41) . lie on the left of j Q . The Cauchy residues theorem yields (x € [0, 1]) 



P (x, s) e st ds = 0. 



This and (pT2"3"]) imply 



lim 



P(x,s)e st ds 



7i 



< lim 



P(x,v + iR) 



>(y+iR)t 



dw = 0. 



Similar arguments hold for the integral along 7 2 . Therefore, by the Cauchy residues theorem, 
integrals along 7 and 7q are equal and the inversion of the Laplace transformation does not 
depend on the choice of sq as well as on the choice of s' Q . 
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The Cauchy residues theorem yields (x £ [0, 1]) 

P (x, s) e st ds = 0, 



where f — 7 U r r (see figure I5T5|1 . with the assumption that all poles, i.e. solutions of (|3.4|l . lie on 
the left of 7 . Let r (R) = y ' R 2 + s 2 ,. Consider 



Ims 



Y'o Yo 

An *S r 



Res 



Ims 




Figure 3.2: Integration contour T 



Figure 3.3: Integration contour f 



P{x,s)e st ds 



tfo(K-R)) 
-4> (r(R)) 



P x, Ji? 2 + s 2 c 



2Jt 



^ R 2 + si cos <f>^± 



where lim </> (r (i?)) = |. Since M ~ y^ as |s| — >• 00, by (I2.9|) and (|3.1|> . we have, as |s| — > oo, 



P(x,s) 
By ([3~24]l . we have 



-(l-i)sM(s) 



i _ „-2isM(s) 



J _ g-2sM(s 



< C, x € [0, 1] , s € C\ (-oo, 0] . (3.24) 



lim 



P(x,s)e s *ds 



•0o('-(-R)) 

< C lim / ,/i? 2 + sgeV **+•<> cos *d</> = 0, 

7?.— s-oo 



since £ < and cos > 0. Therefore, we proved (13.161) . 

By the use of (|3.15[) and (|3.16[) in (I3.14p and by calculating the convolution we obtain (|3.18l) 
and (J3TT5I) . 

In order to prove that uh is a continuous function at t = 0, we will use Lebesgue dominated 
convergence theorem. Let 

,(,,,) := li ( **!?"<■*:*» - ^X ) ■ »C(0,co>,«|0,l|. 



2?ri \ sinh(gM((7e- i7r )) sinh (qM (qe i7r )) 



Then 



1 - c~ qt 
f (q, x) dq 

q 



as t -t 0. 



(3.25) 



By simple calculations we have — - — < Ct if < q < 1 and — - — < 1 — e 9 * if g > 1 . Thus 

f{q,x) l ~ e 9 <Cf(q,x), q>0 
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and since 



as t -> 0, (j3T23t follows. 



In proving the continuity of uh at t = 0, by (|3.14|) . we estimated 

/ P(»,r)dr, i£ [0,1], i > 0, 
Jo 



and actually proved that P is integrable function on any interval [0,1*] , T > 0. Thus, P is locally 
integrable on R. D 

3.1.2 CaseT = T Q H + F 

Condition 3.4. Let F be a locally integrable function, equal to zero for t < 0, smc/i i/iai its Laplace 
transformation exists in C\ (— oo,0] . Assume: 

(i) F is analytic and F ^ in C\ (— oo, 0] ; 

(it) for some a > 1, F (s) ~ |4n-, stC\ (— oo, 0] , as \s\ — > oo; 

(Hi) sF (s) ~ o (1) , 5 e C\ (-oo, 0] , as \s\ -> 0. 



If the boundary condition (12.51) is given by (J3.3I) . then the solution to (|2.3p . (|2.4|l . (|2.5|) . given 
by (|2.8|) in the Laplace domain, reads formally 



u (x, s) = % (#, s) + F (s) P (a;, s) , x € [0, 1] , s e C\ (-oo, 0] , 
and in the time domain it is 

u (x, t) = u H (x, t) + F(t)*P (x, t), x € [0, 1] , t G R. 
The existence of uh is shown in £ )3.1.1[ therefore it remains to show the existence of 

u F {x, t)=F(t)*P {x, t), x € [0, 1] , t € R. 
Let £ > 0. The Cauchy residues theorem yields (x € [0, 1]) 

i u F (x, s) e s 'ds = 27ri ^ fees (u F (x, s) e st , P s { +A + Res (u F (x, s) e'^ps^A , (3.26) 

^ r n=l 

where r = Ti U T 2 U T 3 U T £ U T 4 U T 5 U T 6 U 7 (see figure |3~T|) . Since poles P si ±) of u F are actually 
the poles of P, that are obtained from (I3.4[) and they are simple for n > n , the residues in (|3.26j) 
can be calculated as 



Res 



(u F (x,s) 



e st P s (±) 



- sinh(a;sM(s)) f 

U £[sinh( S M( S ))] 



(3.27) 



„(±) 



The proof that the sum in (|3.26p converges is analog to the one presented in ^3.1.11 
Consider the integral along contour IV It reads 



i'i 



u F (x,s) e st ds 



< 



F(p + iR) P(x,p + iR) 



,(p+ifl)t 



dp. 
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According to (|3.24p and condition 13.41 we have 



lim 

R->oo 



uf (x, s) e st ds 



1 



< C lim 

R ^°° J o ( ^p 2 + i? 2 



« e pt dp = 0. 



The integral along contour I^ reads 
uf (x, s) e st ds 



< 



F (Re 1 ' 



: ,fi(l-a;)c i *M(_Re i '*) 



^2xR.c i *M{B.c i4 ') _ 



3 2Ke 1 <*M(.Re'*) _ I 



^Rt cos 



Rd(j>. 



In order to apply the Lebesgue theorem, we need a > 1 in condition 13.41 (actually it is enough to 
have a > 1, but case a = 1 is already considered). Since M ~ y^ as |s| — > 00 and cos</> < for 
<j) € [f , 7r] , we have 



lim 



/ U_F (X, S) 

Jr 2 



e st d.< 



C lim / i? 1 -« e flcos ^( t+(1 ^ ) \ / *)d0 = O. 



Similar arguments are valid for the integral along the contour IV Thus, 



lim 



The integration along contour T s gives 
uf (x, s) e st ds 



< 



uf (x, s) e st ds 



F (ee i0 ) 



= 0. 



5 (l-x)e I '*M(Ee i *) 



1-e 



-Ixee^M^*) 



Y _ g— 2ec'*M(ee i * 



^ed0, 



and this tends to zero as e — ► 0, according to condition [331 Integrals along parts of contour T3, I\ 
and 7 give 

- , . . sinh(ccgM(oe i7r )) 



hm / uf (x, s) e ds = 

R— >oo /p 

lim / Up (x, s) e st ds - 

R— >oo /p 

lim / iip (x, s) e s 'ds = 27riu^(x,i) 



. , sinh(xqM(oe" 17r )) 

F (qe- 11T ) V • ' e~ qt dq, 

vy ; sinh(gAf (qe-™)) 
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Now, by the Cauchy residues theorem (|3. 261) . u F is determined as 



u F (x,t) = — / [F(qc 



smh. (xqM (qe-™)) 
sinh (qM (qe-™)) 



+ 



- , . . sinh(xM(ge i7r ))\ 
v y sinh(gM(ge 1,r )) / 

J3 [Res (u F e st , p4 +) ) + R es (^e st , p4 _) 



x e [o, l] , t > o, 



(3.28) 



u F (x,t) = 0, aj € [0, 1] , t < 0, 



where the residues are given by (|3.2T|) . Note that up is a locally integrable, real- valued function, 
which can shown similarly as in £13.1.11 

Therefore, in the case when boundary condition takes the form (|3.3|) . the solution to system 
(I23| . ASH), J2HJ reads 

u (x, t) — uh (x, t) + up (x, t) , x € [0, 1] , t > 0, 



where it# and ttj? are given by (J3.18I) and (|3.28[) . respectively. Note that uh and u F are equal to 
zero for t < 0. Again, we have that u is a smooth function for x e [0, 1] , £ > 0. 



3.2 Determination of stress a in case of stress relaxation 



We see that f, given by (|2.13|) . has the branch point at s — and poles at the same points as P. 
Therefore, the poles of T are given as solutions to Q3.4p . Using the Cauchy residues theorem 



(3.29) 



I f (x, s) e st ds = 2iri ^ fees (f (x, s) e st , P s ( +A + Res Of (x, s) e st , P s { ~ 



where contour T is given in figure I3.1[ we obtain T in the following way. Residues in Q3.29J) are 
given by 



Res 



(f(x,s)e st , P s^) 



cosh (xsM (s)) st 

Af( S )£[sinh( S M(s))] e 



(3.30) 



,(± 



where P Sn , n € N, are solutions of (|3.4[) . 
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1 r°° ( cosh (xqM (ge i7r )) 
27ri7 I M (ge i7r ) sinh (qM (qe i7r )) 



Evaluating the integral at the left hand side of (|3.29[) in the same way as in £13.1.11 we obtain 

T(x,t) 

cosh (xqM (ge-™)) \ e 
M (qe~^) sinh (qM (qe- i7r )) J 

+ ^2[Res(f(x, s )e st :P s^) 

+ Rcs(f (x,s)e st , P s^A , x G [0,1], t > 0, 

r(i,t) = o, jc€ [o,i], t<o, 

where the residues are given by (13.301) . The proof is analogue to the one presented in M3.1.H 
Thus, by (J2~T3|) , we have 

a ff {x,t) = ToT(x,t), x£ [0,1], t > 0, (3.31) 



if boundary conditions (|2.5[) are given by (|3.2[) . Note that cr# is a locally integrablc function with 
the jump at t = and smooth for t > 0. Also in case when boundary conditions (|2.5|) are given by 
(O, we have (i€ [0,1]) 

a F {x,t) = a H (x,t) + —(F(t)*T(x,t)), t>0 
a F (x,t) = 0, t<0. 

This is a smooth function for t > 0. Note that an and <jf are real-valued functions, which can 
shown similarly as in ^3.1.11 

4 Numerical examples for displacement and stress in case of 
stress relaxation 

In this section, we give several numerical examples of displacement ur and stress an , given by (J3.18D 
and (|3.31[) . respectively. In figure [4~T1 we show displacements, determined according to (|3.18p . for 
two different positions. Parameters in (|3.18[) are chosen as follows: Yq = 1, a = 0.045, b = 0.5. The 
integration goes to 1000, while the number of residues in the sum is 400. 

In figure 14.21 we show the stresses determined according to (|3.31[) for the same values of pa- 
rameters used for figure |4~T1 In order to emphasize stress relaxation process, we show stresses only 
for t > 1. From figure [4~2l it is evident that the stress tends to a constant value To, at each point 
X G [0, 1] . 

5 Conclusion 

In this work we analyze displacements u and stresses a for a viscoelastic rod of finite length, which 
satisfy a constitutive equation of distributed fractional order (|1.4|) . (|1.5|) . Displacement of a free 
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U H (*.') 





Figure 4.1: Displacements uh (x, t) in stress relaxation experiment, as a function of time at x = 0.25, 
x = 0.75 for te (1,10). 




Figure 4.2: Stresses cr# (a;, t) in stress relaxation experiment, as a function of time at x = 0.25, 
£ = 0.75 for te (1,15). 

end of a rod is assumed to be T = YqH (To > 0) and the displacement uh is obtained in the 
form (13.18[) . Results for displacements and stresses are shown in figures 14.11 and 14.21 Figures show 
oscillatory character of both stresses and displacements. Oscillations are damped and for large 
time displacements show linear dependence on the distance of a particle from a fixed end, while 
stresses are approaching to the limiting value independently of the position of the particle. For 
large time stress relaxation curves tend to curves corresponding to quasistatic analysis (see work by 
Atanackovic 2002a and Drozdov 1998). In figures |4~T1 and l4~2l we show displacements and stresses 
for t > 1. 

Finally, let us comment the choice of parameters in constitutive equation (11.41) , (jl.51) . Our choice 
a < b in (|1.5|) is a result of the requirement that entropy inequality is satisfied for each a G (0, 1) 
(see a paper by Bagley & Torvik 1986). The constitutive equation (|1.4|) . (|1.5[) is of a viscoelastic 
type. A generalization of the problem treated here would include effects of viscoinertial type. In 
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that case (|1.4p would be replaced by 

/ (j) 1 (a) Q D?<T{x,t)da = E f (j> 2 (a) D%£ (x, t) da, x G [0,1] , t > 0. (5.1) 

Jo Jo 

For the analysis of a system (|l.lj) . (|1.3|l . (|5.1|) . one could use the type of analysis presented in work 
by Atanackovic et al. (20096) and Atanackovic et al. (2009c). 

Acknowledgement. This research was supported by Ministry of Science projects 144019 (T.M.A. 
and D.Z.) and 144016 (S.P.). 
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